SCHUR DYNAMICS OF THE SCHUR PROCESSES 



Alexei Borodin 

Abstract. We construct discrete time Markov chains that preserve the class of 
Schur processes on partitions and signatures. 

One application is a simple exact sampling algorithm for q"°'" me -distributed skew 
plane partitions with an arbitrary back wall. Another application is a construction of 
Markov chains on infinite Gelfand-Tsetlin schemes that represent deterministic flows 
on the space of extreme characters of the infinite-dimensional unitary group. 
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Introduction 

The Schur processes were introduced in [OR1] as a class of measures on sequences 
of partitions in order to study large random plane partitions with weights propor- 
tional to q volume j < q < 1. The concept generalized that of the Schur measures 
introduced earlier in [Ok]. The asymptotic techniques of [OR1] were developed 
further in [OR2] to study the asymptotics of large skew plane partitions, see also 
[BMRT] . 

The range of applications of the Schur measures and Schur processes expanded 
quickly; apart from random plane partitions they have been applied to harmonic 
analysis on the infinite symmetric group [Ok] , Szego-type formulas for Toeplitz de- 
terminants [BO] , relative Gromov-Witten theory of C* [OP] , random domino tilings 
of the Aztec diamond [J2] , discrete and continuous polynuclear growth processes in 
1+1 dimensions [PS], [Jl], topological string theory [ORV], and so forth. 

The goal of this paper is to define discrete time Markov chains that map Schur 
processes to themselves, possibly modifying the parameters. We also define Markov 
chains on the two-sided Schur processes introduced below; the principal difference 
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of those from the Schur processes is that they live on sequences of signatures that, 
unlike partitions, may have negative parts. 

The dynamics we construct is also 'Schur like'; for example, an evolution of a 
partition or a signature that represents a fixed slice of the (possibly two-sided) 
Schur process is also a (possibly two-sided) Schur process. 

We present two applications of the construction. 

First, we give an exact sampling algorithm for measures of type q volume on skew 
plane partitions. Other sampling algorithms for such measures are known, see 
[BFP] and references therein. However, it seems that the algorithm we suggest is 
simpler; for skew plane partitions with support fitting in A x B box, the algorithm 
consists in sampling no more that AB(B + l)/2 independent one-dimensional geo- 
metric distributions. A short 'code' for the algorithm can be found in Section 7. 
Exact sampling algorithms for boxed plane partitions based on similar ideas were 
constructed in [BG], [BGR]. 

The second application is a construction of Markov chains on infinite Gclfand- 
Tsetlin schemes that preserve the class of Fourier transforms of the extreme charac- 
ters of the infinite-dimensional unitary group, see Section 4 for details. For similar 
developments on the infinite-dimensional orthogonal group see [BK]. 

A special case of the Markov dynamics that we construct has been studied in 
detail in [BF] . One of the goals of this paper is to provide a more general setup (a 
broad class of initial conditions and a multi-parameter family of update rules) for 
large time asymptotic analysis of the dynamics. 

The construction below is based on a formalism developed in [BF], which in its 
turn was based on an idea from [DF]. However, our exposition is self-contained. 

Acknowledgements. This work was supported in part by the NSF grant DMS- 



1. NONNEGATIVE SPECIALIZATIONS OF THE SCHUR FUNCTIONS 

In what follows we use the notation of [M] . 

Let A be the algebra of symmetric functions. A specialization p of A is an algebra 
homomorphism of A to C; we denote the application of p to / G A as f(p). The 
trivial specialization takes value 1 at the constant function 1 e A and takes value 
at any homogeneous / 6 A of degree > 1. 

For two specializations pi and p2 we define their union p = (pi, p 2 ) as the 
specialization defined on Newton power sums via 



Definition 1. We say that a specialization p of A is nonnegative if it takes non- 
negative values on the Schur functions: s\(p) > for any partition A. 

The classification of all nonnegative specializations is a classical result proved 
independently by Aissen, Edrei, Schoenberg, and Whitney [AESW] (see also [E]) 
and Thoma [T]. It says that a specialization p is nonnegative if and only if the 
generating function of the images of complete homogeneous functions has the form 



0707163. 



Pn{Pl,p2) = Pn(pl) +Pn(P2), 



n > 1. 



(1) 
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for certain nonnegative {ai}, and 7 such that + A) < 00 ■ 

It turns out that nonnegativity of sa(p) for all A is equivalent to nonnegativity 
of the images of the skew Schur functions sx/^p) for all A and p. Hence, via the 
Jacobi-Trudi formula 

sa/p = det[/j Ai -i- w +j]r ij= i, r > max{^(A), %)}, 

the classification of nonnegative specializations is equivalent to that of totally non- 
negative triangular Toeplitz matrices with diagonal entries equal to 1. An excellent 
exposition of deep relations of this classification result to representation theory of 
the infinite symmetric group can be found in Kerov's book [K]. 

For a single a or a single j3 specialization, the values of skew Schur functions are 
easy to compute: 

(2) 

„, x 1 • r , s J a |A| " M > Ai > m > A 2 > M 2>..., 

H{p;u) = - implies s x/ (p) = i 

1 — au 10, otherwise; 

(3) 

H(p; u) = 1 + j3u implies Sx/^p) ■ 



/3l A l"W ; A 3 - p, e {0, 1} for all j > 1, 
0, otherwise. 



We say that a nonnegative specialization p of A is admissible if the generating 
function (1) is holomorphic in a disc D r = {u G C | \u\ < r} with r > 1. In other 
words, p is admissible iff on < r _1 < 1 for alii. 

Since H(p±, P2]u) = H(pi;u)H(p2',u), the union of admissible specializations is 
admissible (unions of nonnegative specializations are also nonnegative). 

For a nonnegative specialization p, denote by Y(p) the set of partitions (or Young 
diagrams) A such that s\(p) > 0. We also call Y(p) the support of p. The set of all 
partitions will be denoted as Y. 

Using the combinatorial formula for the Schur functions [M, Sect. 1.5 (5.12)] and 
the involution u> [M, Sect. 1.2], it is not hard to show that if, for a nonnegative 
specialization p, in (1) 7 = and there are p < 00 nonzero aj's and q < 00 nonzero 
/3j's, then Y(p) consists of the Young diagrams that fit into the T-shaped figure 
with p rows and q columns. Otherwise it is easy to see that Y(p) = Y. 

In particular, if in (1) all /3j's and 7 vanish, and there are p nonzero aj's, then 
Y(p) consists of Young diagrams with no more than p rows. Such a specialization 
consists in assigning values aj to p of the symmetric variables used to define A, and 
0's to all the other symmetric variables. 

We will also need minors of arbitrary (not necessarily triangular) doubly-infinite 
totally nonnegative Toeplitz matrices. The classification of such matrices was ob- 
tained by Edrei in [E], who proved an earlier conjecture of Schocnberg. The result 
is as follows. 

A matrix M — [Mi_j]f?Z_ 00 is totally nonnegative if and only if, after a trans- 
formation of the form M n i-> cR n M n with c > 0, R > 0, the generating function of 
its entries has the form 

+00 

(4) H(M;u):= ^ M n u n 

n— — 00 

= e7 +(„-i )+7 -( u -i-i) TT / ! + #(«-!) l + (3-(u^-l) \ 
l\\l-a+(u-l) l-a7(u-i-l)J 
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for certain nonnegative {a^}, {fif}, and ^ such that { a f + a i + fit + ft 7 ) < 00 
and /3j < 1 for all j. The parametrization of M by {{af}, {fif}i 1^) becomes 
unique if one adds the condition max :( -{/3+} + ma,Xj{(3j} < 1. 

The generating function on the left is understood as the Laurent series of the 
holomorphic function in a neighborhood of the unit circle \u\ — 1 that stands on 
the right. We call the largest annulus of the form {u € C | < r\ < \u\ < r 2 } where 
H(M; u) is holomorphic (the unit circle must be inside the annulus) the analyticity 
annulus of H(M;u). 

Definition 2. We say that a totally nonnegative Tocplitz matrix M is admissible 
if the generating function of its entries is given by (4) (i.e., no multiplication by 
cR n is involved). 

Note that since multiplying Toeplitz matrices corresponds to multiplying the 
generating functions (4), the product of two admissible matrices is admissible. 

It will be convenient for us to use a similar notation for the minors of general 
Toeplitz matrices as in the triangular case (Jacobi-Trudi formula). 

Define signatures of length n as n-tuples A = (Ai > A2 > • • • > A„) of nonin- 
creasing integers. We will also write £(A) = n and |A| = Ai + A 2 + • • • + A„. By 
convention, there is a unique signature of length with = 0. 

For any two signatures A and p of length n and an admissible M we set 



s\/n(M) = dot [M Xi _i_ H+j \ 



For totally nonnegative M with only one or (3^ parameter nonzero (and all 
other parameters being zero), one obtains formulas analogous to (2), (3): 

1 1 / a \ ±1^1 

(5) g fo u > = i-« (u ±i-i) implies = (TT^ V T+W 

if ±Xj =F > f° r all 1 < j < n, and otherwise; 

/ a \±|A|tImI 

(6) H(p;u) = l + f3(u ±1 -1) implies s A/ » = (1 - p) n (t^J 

if ±Xj =F pj € {0, 1} for all 1 < j < n, and otherwise. 

Also, mimicking the property of the Schur functions, for a constant c e C, a 
signature v of length n + 1, and a signature A of length n, we set 

(?) ^ (c) . = f cM-M, A„ +1 <Mn<A n <---<A 2 <Mi<Ai, 

\ 0, otherwise, 

with the convention that 0° = 1. 



2. The Schur processes 

Pick a natural number TV and admissible specializations p$ , . . . , p^_ 1 , p{ , . . . , p^ 
of A. For any sequences A = (A^, . . . , A^) and p = (p^\ . . . , p( N ~^) of parti- 
tions satisfying 

(8) c AW d m « c A< 2 > d m ( 2 ) c ■ ■ o p^-V c A< w > d 
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define their weight as 
(9) 

W(A, p) := s A( i) (p+) s A (i )/m (i) (p x )s A (2 )/aj( i) [pf) ■ ■ • s A («) /ai (jv-i) (pj i _ 1 ) s xm (p N ). 

There is one Schur function factor for any two neighboring partitions in (8). 
The fact that all the specializations are nonncgative implies that all the weights 
are nonncgative. The admissibility of p's implies that 

(10) Z(p+,...,p+_ i;P i,...,p N ):=J2 W ( X ^)= II H(p+;pJ)<<x>, 

X,H 0<i<j<N 

where H(p 1 ;p 2 ) = Easy s a(pi)sa(P2) = cxp (^2 n > 1 p n (pi)Pn(P2)/nj, andp„'s are 
the Newton power sums. Indeed, this follows from the repeated use of identities, 
cf. [M, 1(5.9) and Ex. 1.5.26(1)], 

(11) X] S */v(Pl) S K/i>(p2) = H (P1:P2) X s v/t(P2)s0/ t (pi), 

kGY tGY 

(12) S k/"(Pi) s v/t(P2) = S K / T (p 1 ,p 2 ), 

and from the fact that for an admissible specialization p with H(p; u) holomorphic 
in a disc of radius r, we have p n (p) = 0(r~ n ). 

The same argument shows that the partition function (10) is finite under the 
weaker assumption of finiteness of all H(p+; pj) for < i < j < N. 

Definition 3. The Schur process S(pq , . . . , p%_i, P\ , ■ ■ ■ , p^) is the probability 
distribution on sequences (A, p) as in (8) with 

c + + —\(\ \ W(A,p) 

S(Po i ■■•■> Pn-vPi > • • • . Pat)(A, P) - 



Z{Pa,---,P%_ 1 ]P 1 ,---,P N ) 

The Schur process with N — 1 is called the Schur measure. 

Using (11)-(12) it is not difficult to show that a projection of the Schur process 
to any subsequences of (A, p) is a also a Schur process. In particular, the projection 
of S(p^,...,pj { _ 1 ;pi,...,p N ) to AW is the Schur measure S(p^ )J ._ 1] ;py |JV] ), and 

its projection to p^ is a slightly different Schur measure S(p^ fc _ 1 p p^ +1 N ^)- Here 

we used the notation p^ b , to denote the union of specializations p^j, m = a, . . . ,b. 

We now aim at defining a Schur like process for signatures. 

Pick a natural number AT, real numbers oi,...,ajv > 0, nonnegative integers 
c(l), . . . , c(N), and c(l) H + c(iV) admissible Toeplitz matrices 

M = {M (fc < z) | 1 < k < N, 1 < I < c(k)}. 

If all c(k) are zero then M is empty. 

We will also need a totally nonnegative matrix of size ZxiV, denote it as 

* = [**j]»ez,-i>j>-jv 
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For any sequences 

(13) A« = (A* 1 - ), A' 1 - 1 ', . . . , A^ 1 ))), . . . , \W = ( A W), A^ 1 ', . . . , A^ C W)) 
of signatures of lengths £(A^ fe '*^) = k, define their (nonnegative) weight as 

(14) W(A (1) ,...,A (Ar) ) :=det < 



-1,-3 



N 



AT 



c(fc) 



0) /^(*-l,c(k-l)) 



fe=l 



with A(°- c (°» = 0. 

We assume that the generating functions 

(15) *,-(«):= >J *„,_,- ^ 



n— — oo 



are holomorphic in an open set containing the unit circle. As we will see in 
Section 10, if for any j < N, aj lies in the common analyticity annulus for 
{^(MWlju- 1 )}^, then the partition function of weights (14) is 

finite and it has the form 

(16) Z(a u ...,a N ;M;*):= ]T W(A« , . . . , A< s ') 

At 1 ),..., A< N > 

AT 

c(fc) 



det J ^j(ai) 



det 



JV 



i^i n u H ( Mik ' i w)- 



l<j<k<N 1=1 



In the important special case when the matrix <]/ is actually Toeplitz, ^i.-j = 
r/'i+jj (16) simplifies: 

N c(fc) 

(17) Z(oi,...,a JV ;M;*) = JJ^(oi) JJ IT (M^aT 1 ) , 

i=l l<j<k<N 1=1 

where V(«) = E n ez^« u ™- 

Definition 4. The two-sided Schur process T(ai, . . . , a at; M; ^) is the probability 
distribution on sequences (A^, . . . , A^) as in (13) with 

W(A«,...,AW) 



T(ai,...,a JV ;M;*)(A< 1 >,...,A< JV >) - 



Z(ai 



Remark 5. If in the Schur process of Definition 3 each of the specializations is 
a one- variable specialization with H(pf;u) = (1 — aj + iti) _1 , j = 0, . . . , N — 1, then 
the Schur process can be viewed as a special case of the two-sided Schur process 
with c(l) = • • • = c(N — 1) = 1, c(N) = 0, and identification 

a0')=A^°\ j = l,...,JV, /i^) =A( J ' >1 >, j = l,...,JV-l, 

^; U )=i/(M^ 1 ); U - 1 ), fc = l,...,AT-l; fffe;u) = ^). 

The corresponding two-sided Schur process lives on signatures with nonnegative 
parts that can also be viewed as partitions. 

Observe that under this identification the formulas (10) and (17) coincide. 



SCHUR DYNAMICS OF THE SCHUR PROCESSES 



7 



3. Example 1. Measures q volume on skew plane partitions 

Fix two natural numbers A and B. For a Young diagram 7r C B A 1 set 7f = B A /it. 

A (skew) plane partition II with support fr is a filling of all boxes of n by 
nonnegative integers IT j (we assume that 11,-^ is located in the ith row and jth 
column of B A ) such that EEj j > n^j+i and Il^j > Ilj+i ^ for all values of 

The volume of the plane partition II is defined as 

voi(n) = J2n iJ . 

i,j 

The goal of the section is to explain that the measure on plane partitions with 
given support n and weights proportional to q vo "'\ < q < 1, is a Schur process. 
This fact has been observed and used in [OR1], [OR2], [BMRT] . 

The Schur process will be such that for any two neighboring specializations 
Pk'Pk a ^ l eas t one is trivial. This implies that each fi^ coincides either with A^ 
or with Thus, we can restrict our attention to A^'s only. 

For a plane partition II, we set (1 < k < A + B + 1) 

AW(n) = {n^+^A-i I + k - A - 1) e 7f}. 

Note that A^ 1 ) = A( A+S+1 ) = 0. 

We need one more piece of notation. Define 



£(tt) = {A + m - i + 1 | i = 1, . . . , A}. 

This is an A-point subset in {1, 2, . . . , A + B}, and all such subsets are in bijection 
with the partitions n contained in the box B A . The elements of C(tt) mark the 
"up-steps" in the boundary of tt (=back wall of IT). 




The figure above shows a plane partition IT and its plot with 
A = 4, B = 3, tt= (2,1,1,0), 

A (2) = (4)) A (3) = 3) A (4) = (5) 1}) A (5) = (1Q) 2)) A (6) = (6)) A (7) = (8)) 
A+B 

voi(n) = I aW I = 39 > £ W = i 1 ' 3 - 4 > 6 1- 

i=2 
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Proposition 6. Let it be a partition contained in the box B A . The measure on the 
plane partitions II with support n and weights proportional to q vol ( n \ is the Schur 
process with N = A + B + 1 and nonnegative specializations {pi}, {pj} defined by 

H{p+:u) = H(p N ;u) = 1, 

Note that not all specializations are admissible, but the weaker assumption of 
fmiteness of H(pf; pj) for < i < j < N guarantees that the partition function is 
finite. 

Proof. Observe that the set of all plane partitions supported by 7f, as well as the 
support of the Schur process from the statement of the proposition, consists of 
sequences (A«, A< 2 \ . . . , \W) with 

A (1) = AW = 0, 
\U) \U+V if j g C(X), X U) >■ X u+1) if j i C(\), 

where we write p, ~< v or v >- p iff z^i > pi > v 2 > P-2 > ■ ■ ■ ■ 

On the other hand, (2) implies that the weight of (A (1) , X (2 \ . . . , A (JV) ) with 
respect to the Schur process from the hypothesis is equal to q raised to the power 

A+B 

I a0) |( _ 0' - 1 ) 1 j-ie£(7 r ) - (j - l)l;7-ig£W + .?'lje£(7r) + .?'l^£( 7r )), 

where the four terms are the contributions of pj_ 1 , pj_-^, pj , pj , respectively. 
Clearly, the sum is equal to J2f=^ |A^'| = vol(U). □ 

Remark 7. A similar statement holds for any measure on plane partitions with 
weights proportional to n?]^ with possibly different positive parameters qj, as 
long as the partition function is finite. The proof is very similar. 

4. Example 2. Path measures for extreme characters of U(oo) 

Let U(N) denote the group of N x N unitary matrices. It is a classical result 
that the irreducible representations of U(N) can be paramterized by signatures 
A = (Ai > ... > Xn) of length N also called highest weights. Thus, there is a 
natural bijection A o \ X between signatures of length N and the conventional 
irreducible characters (=traces of irreducible representations) of U(N). 

For each N, embed U(N) in U(N + 1) as the subgroup fixing the (N + l)st basis 
vector. Equivalently, each U £ U(N) can be thought of as an (N + 1) x (N + 1) 
matrix by setting C/^jv+i = Un+ij = for 1 < i,j < N and ?7jv+i,jv+i = 1- The 
union \J^ =1 U(N) is denoted U(oo) and called the infinite- dimensional unitary 
group. 

A character of U(oo) is a positive definite function \ : U(oo) — > C which is 
constant on conjugacy classes and normalized by x(l) = 1. We further assume 
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that x i s continuous on each U(N) C U(oo). The set of all characters of U(oo) is 
convex, and the extreme points of this set are called extreme characters. 

Remarkably, the extreme characters of U(oo) are in one-to-one correspondence 
with admissible Toeplitz matrices M from Definition 2, see [Vo], [VK], [00]. The 
values of the character \ M corresponding to M are given by 

X M (U)= J] H(M;u), 

wGSpcctrum(C7) 

where H(M; u) is given in (4). 

Let GTjv be the set of all signatures of length N; set GT = IJ^ GTjy. Turn GT 
into a graph by drawing an edge between signatures A £ GTjv and fi £ GTjy+i if 
A and [i satisfy the branching relation A -< fi, where A -< fi means that /ii < Ai < 
fJ-2 < A2 < . . . < Xn < fJ-N+i- GT is known as the Gelfand-Tsetlin graph. 

A path in GT, or an infinite Gelfand-Tsetlin scheme, is an infinite sequence 
t = {ti,t2, . . • ) such that ti £ GTi and ti -< ti+i- Let T be the set of all such paths. 

One can also look at finite paths, or finite Gelfand-Tsetlin schemes, which are 
sequences r — {t\,T2, ■ ■ ■ ,tn) such that Ti £ GT^ and t\ < t-i -<...-< t^. Denote 
the set of all paths of length N by 7jv- 





The figure above depicts a Gelfand-Tsetlin scheme r £ Ta and its plot with 

n = (3), 7a = (4,-1), r 3 = (5,0,-5), = (5, 1,-2,-7). 

Each character x of U (00) defines a probability measure Ph on GTtv: Restricting 
the character to U(N), we have 



X 



U(N) 



X 



X A (1jv) 



For each finite path t £ Tn , let C r C T be the set 

C T = {t £ T : (t u t 2 ,... ,t N )=r}. 

A character x of U(oo) also defines a probability measure P x on T (with a 
suitably defined Borel structure), which can be uniquely specified by setting 



P X {C T ) 



i$(A) 



X A (1 



N) 



where r is an arbitrary finite path ending at A, see [01, Section 10] for details. Note 
that we assign the same weight to all finite paths with the same end. 

We use the same formula to define a probability measure P% N ^ on 7jv j which is 
just the projection of P x from T to Tn- 



10 



ALEXEI BORODIN 



Proposition 8. For any admissible Toeplitz matrix M as in Definition 2, the 

M 

measure P^ N ^ on Tn coincides with the two-sided Schur process of Definition 4 
with 

oi = ••• = a N = 1, c(l) = • • • = c(N) = 0, * = M, 
and with sequences (A' 1 ' ), . . . , A^' )) viewed as elements oJTn- 
Proof. Directly follows from (7) and Lemma 6.5 of [01]. □ 

5. Markov chains on the Schur processes 
Let us introduce some notation. 

For two nonnegative specializations pi,p 2 of A such that H(p\;p 2 ) < oo, and 
A, u E Y, set 

Ppi,P2 

(A, \i tt v) = const ■ s v /x{pi)s v/fl {pi), v e Y, 
where we assume that 

(18) {v e Y I s v/x ( Pl )s v/fl (p 2 ) > 0} ^ 0, 

and the constant prefactor is chosen so that we obtain a probability measure in v. 

j2Ppup 2 (^»n u) = i. 

Given (18), the existence of such constant follows from (11). 
Similarly, dropping the assumption H{p\;p 2 ) < oo, we define 

PpuPzi^Vlt v) = const ■ s X / v (pi)s v/lJi {p 2 ), 

Ppi ,P2 

(A, /j,tiv) = const ■ s v /x{pi)s lx / v {p2), 

Ppi ,P2 

(A, ^ LJ. v) = const ■ s X / v (pi)s IJ ,/ l/ (p2), 

where in all three cases we assume that the set of v giving nonzero values on the 
right-hand side is nonempty (it is finite in all three cases), and we choose constants 
so that we obtain probability distributions in v 6 Y. 

If both pi and p 2 are single-a or single-/? specializations, relations (2), (3) show 
that all four distributions P Pl yP2 are products of geometric distributions conditioned 
to stay in segments and Bernoulli measures. 

Example 8. Assume that H(p\;u) = (1 — au)^ 1 , H(p 2 ;u) = (1 — bu)^ 1 . Denote 
by G|j >n , m < n, the probability distribution on the set {m, m + 1, . . . , n) given by 

Then 

PpuPii^it 1 tt V) = Gm ax {A 1)Ml },+oo( I/ l) II Cmax{A 3 , W },min{A J _ 1 ^ 3 _ 1 }( iy i)i 

P P1 P2 (\,filtv) = G b/a , . (vi)T\G h,a /x , . n ,(i/A 

P1,P2\ it t y max{A 2 },Ai v 'XX max{A J + i },min{A 3 } v -' y 

i>2 
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where in the first case we need to additionally assume that ab < 1 (equivalently, 
H(p 1 ;p 2 ) < oo). 

Further, assume that H(p 3 ; u) = (1 + cu). Denote by B^ n , n £ {m, m + 1}, the 
probability distribution on {m, m + 1} given by 

1, k = m, . . . | YTc' k = m ' 



, fc = m + 1. 



1 + c 



Then 

Pp3,pA\ P-W u ) = - B max{Ai, A ii},Ai+l(^l) \\ A 3 , w } ,min{ A 3 _ i , Aj + 1, W - 1 } ( U 3 ) ' 

J">2 

-Pp3,p 2 (^! P 4-t f) = ^maxfAi-l^M.Ai^ 1 ) II ^m'axj Aj -1, X J + 1 , H } ,min{ A 3 i } ("j ) ' 

i>2 



Let (/9q , . . . , pjf_i, Pi, ■ ■ ■ , Pn) be nonnegative specializations of A defining a 
Schur process as in Definition 3. Let n be another nonnegative specialization of A 
such that H(tt, p+) < oo for all < j < N. 

Let X be the set of pairs of sequences (A, p) as in (8) with 

S A(D (Po ) S A(D/ M <i) (Pl)sx(2) /f _ l (i) (pt) ■ ■ ■ SxW/^N-i) O^-i) > 

The product above is the same as in (9) without the last factor. Thus, the support 
of S(pq~, . . . , plf_i, Pi , ■ • • , Pn) is contained in X. 

Define a matrix CpJ. with rows and columns parameterized by elements of X via 

(19) <pt((A, p), (A, p)) = P,+ „ (0, A« n A«) 

x U P P7 J~ x{j) ^ U) ^A W) )i > ,+ w (A W) ,A«+ 1 ) ttA^ 1 '). 

In other words, starting from (\,p), one first finds A^ using A^, then /iW 
using A^ and p^\ then -V 2 ) using p^ and A^ 2 \ and so on. One could say that 
we perform sequential update. 

Note that some of the entries of might remain undefined if one of the con- 
ditions of type (18) is not satisfied. Part of the theorem below is that this never 
happens. 

Theorem 10. In the above assumptions, the matrix ty\ is well-defined and it is 
stochastic. Moreover, 

s (Pct, • • -,Pn-i,Pi, ■■■,Pn)VI = s (Pct, • • -,Pn-i;Pi, ■ ■ -Jpn), 

where "p" = (p^,7r). In other words, changes the last specialization of the 
Schur process by adding n to it. 

The proof of Theorem 10 will be given in Section 9. 
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Matrices <p£ describe a certain growth process. In a similar fashion, one obtains 
a process of decay. Let us describe it. 

Let a be a nonnegative specialization of A that 'divides' , that is, there exists 
a nonnegative specialization b p$ such that p^ = ( b p^,a). For example, a may 
coincide with pj; in that case b p$ is trivial. 

Let y be the set of pairs of sequences (A, p) as in (8) with 

s a(d Cpo) s a(d/ m (i) (prK^/pd) 0i~) ■ ■ ■ sxm/^N-i) > 0. 

Note that if a = p£ then A' 1 ) and p^ must be empty in order for (A, p) to lie 
in y. 

Define a matrix *p£ with rows parameterized by X and columns parameterized 
by y via 

q3i((A, M ),(A, / i)) = P po+!(T (0,A( 1 ) UA«) 

JV-1 

Notice that the only difference of this definition and that of CpJ above, is switching 
7r and a and changing the second arrows from t to 4- 

Theorem 11. Ira i/ie a&cwe assumptions, the matrix ^ is well-defined and it is 
stochastic. Moreover, 

s (Po> ■ ■ -iPn-vPii ■ ■ ■ ■> Pn)^ = s ( b Po,---,P~N-i'-,Pi, ■■■,Pn), 

where ( b pQ,a) = Pq ■ In other words, changes the first specialization of the 
Schur process by removing a from it. 

The proof of Theorem 11 will also be given in Section 9. 

Remark 12. Both Theorems 10 and 11 can be generalized as follows. Assume we 
have an arbitrary sequence of Markov steps of types ffi and ^ applied to an initial 
Schur process, and let us denote by (X(t),p(t)) the result of the application of t 
first members of the sequence. One can show that any finite sequence of random 
partitions of the form 

(\W(t hl ), A«(i 1>2 ), . . . (*;,!), m (1) (*'i, 2 ), 

■ ■ ■ , ^ (JV - 1) (4-i,i).M (JV " 1) (4-i 1 2), ■ ■ • , A (N) (^,i), AWfe), . . . ) 
forms a Schur process with an explicitly known specializations as long as 

tlA > £1,2 > • • > £'1,1 > t' 12 > • • > ijv-i,i ^ *JV-1,2 > ' " ' > £jv,i > ijv,2 > • ■ • , 
cf. the last sentence of Section 8. 
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6. Markov chains on the two-sided Schur processes 

We now aim at formulating (and later proving) a statement for the two-sided 
Schur processes that is analogous to Theorem 10. 

For two admissible matrices Mi and M 2 ('admissible' is explained in Definition 
2), and two signatures A and u of length n > 1, we define a probability distribution 
on GT„ (=the set of all signatures of length n) via 

f , M 1 ,M 2 (A, fi || v) = const ■ s v /x(M 1 )s v / li (M 2 ), v G GT„. 

For an admissible matrix M and a positive number a in the annulus of analyticity 
of H(M; u), and for two signatures A <E GT„_i and fi G <GT„, we define a probability 
distribution on GT n via 

Pa,Ai(\ A* II v) = const ■ s„/x(a)s„/^(M), v G GT„. 

In both definitions, we suppose that the set of z/s giving nonzero contributions 
to the right-hand sides is nonempty. Then our assumptions imply the existence of 
the normalizing constants. 

Similarly to the one-sided Schur process, if M\ and M 2 are both single-a ± or 
single- {3^- matrices, then Pmi,m 2 splits into a product of geometric/Bernoulli ran- 
dom variables, cf. (5)-(6) and Example 8. For P a ,M the same holds if M is a 
singlc-a^ or single- ^ matrix. 

Consider the two-sided Schur process of Definition 4, and let 

X= |(A«,...,AW) G (<GT 1 ) c ( 1 >+ 1 x •■• x {GT N ) c{N)+1 | 

N c(fc) ~j 

]] (sj(M)/J(n,.(H» J] S A (M) / A(M-l) (M ,W) ) >0 ) (' 

fe=l J=l J 

where A(°< c (°» = 0, cf. (14). Clearly, supp T(a u . . . , a N ; M; *) C A". 

Let Q be an additional admissible matrix such that all the parameters Oj lie in 
the analyticity annulus of H(Q;u). Define a matrix *}3q with rows and columns 
parameterized by X via 

*P Q ((A( 1 ),...,AW),(/rW,...,/rW)) = 

AT / c(fc) 

n ^.o (a^- 1 ^- 1 )),^ ) hm (m) ) n p M(M) ;Q (aw-",^') ii^)) 
k=i \ i=i 

The structure of is such that to compute its row indexed by (A^, . . . , 
one first finds /Z 1 ' ) using A^ 1 ' -*, then using A^ 1 ' 1 ) and /j,( 1,0 \ then using 

A' 1,2 ) and /t/ 1 ' 1 ), and so on. 

Theorem 13. In the above assumptions, the matrix Cpg is well-defined and it is 
stochastic. Moreover, 

T(ai, . . . , a N ; M; tf) <£ Q = T(ai, . ..,a N ; M; Qtf ). 



The proof of Theorem 13 will be given in Section 10. 



14 



ALEXEI BORODIN 



Remark 14. Similarly to Remark 12, a more general statement can be proved. 
Assume we have an arbitrary sequence of matrices CPq applied to a two-sided Schur 
process T(ai, . . . , ajy; M; Denote by (X^\t), . . . \( N \t)) the random sequence 
obtained after the application of t first matrices. Then any sequence {^ k,l \tk,i)} 
forms (a marginal of) an explicitly describable two-sided Schur process as long as 
(ki,h) < (k2,h) lexicographically implies t klill > tk 2 ,i 2 - 

Remark 15. The matrices *Pq are similar to the growth process defined by of 
the previous section. One could also define a 'decay process' for the two-sided Schur 
processes that would be similar to ; the application of the corresponding matrix 

to T(oi, . . . , a at; M; ^) would reduce N by 1 and remove a\ and {M^Y^l from 
the set of parameters. 

Remark 16. In the setting of Remark 5, one easily shows that and *}3q coincide 
ii H(n;u) = H(Q;u). 

7. Exact sampling algorithms 

Let us start with (one-sided) Schur processes. Theorem 10 yields an exact sam- 
pling algorithm that is inductive in N. 

As the base one can take the empty sequence and N — 0. Let us explain the 
induction step. Assume we already know how to sample from the Schur process 
V n -i = S(p£,...,pj T _ 2 ;p7;,...,px_ 1 ). 

Consider the process V n = S(pJ, . . . , p%_ 2 , Aw-i! Pi, ■ ■ ■ > Pjv-i' ®)-> where 
is the trivial specialization. The definition of the Schur process implies that for 
this process A^ = p^ N ^^ = with probability 1, and the distribution of the 
remaining partitions (A^ 1 ), pP^, ■ ■ ■ , p^ N ~ 2 \ A^ -1 )) is the same as for V n -\ that 
we already know how to sample from by the induction hypothesis. 

In order to obtain a sample of V n = S(pJ, . . . , p^_ 2 , Pn-i, Pi, ■ ■ ■ , Pn-i, Pn) 
we apply the stochastic matrix ty\ with 7r = p^ to V n , cf. Theorem 10. The 
application of this matrix requires sequential update from A*- 1 -* up, cf. (19). 

We thus see that if each of (p^ , . . . , p^_ 2 , pti-v P\ ■>•••■> Pn -it Pn) * s a singlc-a 
or a single-/? specializations (or trivial), then exact sampling is reduced to sampling 
a finite number of independent geometric/Bernoulli random variables. Noting that 
in the algorithm for the iVth step one docs not have to use a single ty\ with 
7r = p^, but can instead use a sequence of ty\ with p^ = (m, n2, ■ ■ ■ ), we see that 
the a similar reduction holds for the Schur processes with all specializations having 
finitely many nonzero a's and j3's (and 7 = 0). 

For the measures q volume n skew plane partitions considered in Section 3, the 
algorithm can be implemented as follows (we use Section 3 and Example 8 below). 

Initiate by assigning A (1) = • • • = \( A + B ) = 0. 
For k running from 2 to (A + B) 
If k <£ £(tt) then 

For 1 running from 1 to (k — 1) 

If I e C(tt) then A^ +1 ) := v with v distributed as 

r\(0 i v^ 1 ) 11 r\(0 • r\(0 vC+ilirj) 

max{Ai ',X\ T '}, + oo - >2 max{ X\> , \\ T },min{ AV 1 , AS J 1 '} 

If I £ C(tt) then A^ +1 ) := v with v distributed as 
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G + s(l + l)( V l) II G r,(I) ,(1 + 1)-, ■ ,,(!) .(l + lJiC^j) 

max{A^',Ai '},A^ max{A^ 1: A) },mm{A) , A) J ± } 

End of ^-cycle 
End of fc-cycle 

At the end of each fc-step we see an exact sample of the measure q volume on 
plane partitions with a smaller support. The number of nontrivial one-dimensional 
samples needed to go through the k-step with k £ C(n) is the number of boxes 
in this support. It is not difficult to see that this number is at most A for the 
smallest k ^ £(71"), it is at most 2 A for the next one and so on, so that the total 
number of one-dimensional samples needed is at most AB(B + l)/2. The maximum 
is achieved at C(tt) = {1, . . . , A}, i.e. when the plane partitions are supported by 
the full Ax B box. 




The above figures show a sample for a specific back wall profile, and an average 
over ten samples with the same back wall. A limit shape and its cusp are clearly 
visible, cf. [OR2]. 

Finally, note that a very similar algorithm would sample skew plane partitions 

with weights of the form Y[ ' • 

Let us now discuss the two-sided Schur process. First, let us restrict ourselves 
to the case when is Toeplitz. Then if all H(M^ k ' 1 ^ ; and ip(u) are analytic 
in a disc of radius > 1 (not just in an annulus containing the unit circle), then 
the two-sided Schur process lives on signatures with nonnegative coordinates and 
it constitutes a special case of the (one-sided) Schur process, cf. Remark 5. Con- 
sequently, if all M^ k ^ and \P are admissible matrices with M^ k ' 1 ' having finitely 
many a~ and j3~ nonzero parameters (all others are zero), and ^ having finitely 
many a + and /3 + nonzero parameters, the inductive algorithm for the Schur pro- 
cess described above reduces sampling to a finite number of independent samples 
of geometric/Bernoulli random variables. 

On the other hand, Theorem 13 allows us to add finitely many a ± and ^ 
parameters to ^ by sampling from independent geometric/Bernoulli distributions. 
Hence, we can relax the assumption on ^ in the previous paragraph by requiring 
that it has finitely many and (3^ parameters. 
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The figure above shows a sample of the path measure and the average over ten 
samples for the extreme character of U(oo) with 

= = -jj, # = ..-# = 5, a- = --- = a- = ^, 

and all other parameters being zero, cf. Section 4. The first order asymptotic 
behavior of such measures as the path length goes to infinity and parameters remain 
fixed is known, see [00]. 

8. A GENERAL CONSTRUCTION OF MULTIVARIATE MARKOV CHAINS 

The general construction of this section will be used in Sections 9 and 10 to 
prove Theorems 10, 11, and 13. 

Let (Si, . . . , S n ) and («Si, . . . ,«S n ) be two n-tuples of discrete countable sets, 
Pi , . . . , P n be stochastic matrices defining Markov chains Sj — > Sa . Also let Af , . . . , 
A"_ x and Af , . . . , A"_ x be stochastic links between these sets: 

P k ; S k x S k -> [0, 1], J~] Pk(x,v) = 1> xeSk, fc=l,...,?i; 
yes k 

A^_ x : S k x S k -i -»• [0, 1], Al_i(x,y) = 1, iG5 fc , fe = 2,...,n; 

yeSk-i 

A*_i : <5fe x <Sfe_i ->■ [0, 1], £ A^_ 1 (^ 2/ ) = l ! i£S t) fc = 2,...,n. 

Assume that these matrices satisfy the commutation relations 
(20) Ati:=AtiiW=PfeAti, fc = 2,...,n. 

We will define a multivariate Markov chain PW between the state spaces 

n 

S {n) = {(x u ...,x n ) €Si x ••• xS n I [] A^^Xfe-O ^0} 

fc=2 
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and 

n 

S (n) = {(x u ...,x n ) e5i x ••• xS n I JjA^^Sfc.Sfc-O^o}. 



fc=2 



The transition probabilities for the Markov chain P^™) are defined as (we use the 
notation X n = (xi, . . .,x n ), Y n = (yi, . . . ,y n )) 

(2i) p( '(x rt) y„) = p l( , 1 , yi) n aO^M 



if Ilfe=2 ^fe-i( a; fc, 2/fc-i) > 0, and otherwise. 
One way to think of P( n ) is as follows. 

Starting from X = (x\, . . . ,x n ), we first choose y\ according to the transition 

matrix P 1 (x 1 ,y 1 ), then choose y 2 using P2 ^^(if y^) 2 '^ , which is the conditional 

distribution of the middle point in the successive application of P 2 and Af provided 
that we start at x 2 and finish at yi, after that we choose yz using the conditional 
distribution of the middle point in the successive application of P3 and A 2 provided 
that we start at x 3 and finish at y 2 , and so on. Thus, one could say that Y is 
obtained from X by the sequential update. 

Proposition 17. Let m n be a probability measure on S n . Let rnS n > be a probability 
measure on defined by 

m {n) {X n ) = rn n {x n )K^_ 1 {x n ,x n -i)---\\{x 2 ,x{), X n = (x u . . . ,x n ) e S {n) . 
Set rh n — m n P n and 

m {n) {X n ) = TO„(x„)A"_ 1 (x„,x n _i)--- k\{x 2 ,x{), X n = (x u . . . , x n ) e S {n) . 
Then m(")p(™) = mW. 

Proof. The argument is straightforward. Indeed, 



m WpW(y n )= m n (x n )k n n _ 1 {x n ,x n - 1 )---k\(x 2 ,x l ) 

A Pc(x fc ,2/fc)A^_ 1 (y fe ,y fc _i) 



Extending the sum to £1 6 S\ adds to the right-hand side. Then we can use 
relation (20) to compute the sum over x\, removing A^(x 2 ,xi), P\(x\,y\) and 
A-i(x 2 ,yi) from the expression. Similarly, we sum consecutively over x 2 ,...,x n , 
and this gives the needed result. □ 

Proposition 17 will be used to prove Theorems 10, 11, and 13. A more general 
[BF, Proposition 2.7] is needed to prove the statements mentioned in Remarks 12 
and 14. 
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9. Application to the Schur processes 
In this section we prove Theorems 10 and 11. 

Let us start by putting the Schur process S(pg , . . . , p%_ x ; Pi , ■ ■ ■ , p^) of Defini- 
tion 3 into the framework of the previous section. 
We need some general definitions. 
Let y, z, t be nonnegative specializations of A. Set 

where Y(p) = {k E Y s K (p) > 0}, and for the first definition we assume that 

H {y\ z ) = E KeY s k(2/K( 2 ) < °o- 

Relations (11) and (12) imply that the matrices 

p*(y;z) = [pl^y;z)] x ^ eY{y) and p±{y;t) = [p{ v (y; t)] XeYM ^ Y{y) 
are stochastic: _ 

E pUv> z )= E pL(»;*) = 1 - 

It is immediate to see that p^ and p^ act well on the Schur measures: 

(22) S{x; y)p [ {y- 1 z) = S{x, z; y), S{x; y, i)p% t) = S(x; y). 

Observe that S(pi;p 2 ) = S(p 2 ;pi), so the parameters of the Schur measures in 
these relations can also be permuted. 

Proposition 18. Let y, z, z\, z 2 , ti, ti be nonnegative specializations of A. Then 
we have the commutativity relations 

pt(y; z 1 )p i (y- 1 z 2 ) = p t (y; z 2 )p t (y; zi), 
P l (y, t 2 ;t- L )p l (y; t 2 ) = p l (y, t 1 ;t 2 )p i (y; h), 
P t (y, t; z)p l (y; t) = p i {y; i)p r (y; z), 

where for the first relation we assume H(y; z\,z 2 ) < oo, and for the third relation 
we assume H(y, t; z) < oo. 

Proof. The arguments for all three identities are similar; we only give the proof of 
the third one which is in a way the hardest. We have 

1 s„(y) m—^ H(t;z) s u (y) ^ 

= h^.m^'^MV = H{yMz) SA(Sit) g»v.('K/.W 
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where along the way we extended the summation in p from Y(y, t) to Y because 
s v (y)s^/ v {t) > implies s ll (y, t) > by (12); we used (11) to switch from p to k, and 
finally we restricted the summation in k from Y to Y(y) because s v {y)s v / K {z) > 
implies k C v and s K (y) > 0. □ 

We are now ready to return to the Schur process S(pq , . . . , p%_i, P\ , ■ ■ ■ , p~^)- 
Set n = 2N - 1 and 

S 2j -i=Y(p+ t ._ 1] ), j = l,...,N; 
S 2k =Y{p+ k _ 1] ), fc = l,...,iV-l. 

Since and p^ are distributed according to the Schur measures S(p^ ; jv]) 
and S(pi P[k+i at] ) respectively, the projections of the support of the Schur 
process to these coordinates lie inside <S2j-i an d S2k, respectively. 
Dehne the stochastic links by 

A| +1 =^(p [ f 0j _ 1] ;p+), j 1 V 1: 

^j-i=P*{pV,j-iYP7)> j = l,-..,N-l. 
One immediately verifies the formula 
(23) S(p+, . . . , p+_ 1 ;pi,. . . , p N )(X, p) 

= S (p+^j ; p-) (A W ) nfA^ 1 (A^ , „ (fe) ) Alt, (M (fc) ,A«)), 
fe=i ^ ' 

cf. the dchnition of in Proposition 17. 

Proof of Theorem 10. We apply Proposition 17. Set <Sj = Sj for j = l,...,n, 
Aj_ 1 = A^._ 1 for j = 2, . . . , n, and also 

W„ = S(pp JV _ 1 pP JV ), 

f2j-i =j) T (p^._ 1] ;7r), j = l,...,iV; 

=p t (p^_ 1] ;7t), j = i,...,JV-i. 

The commutation relations (20) follow from Proposition 18, and the matrix of 
transition probabilities P^ from (21) is easily seen to coincide with The 
claim now follows from (23), Proposition 17, and the relation (cf. (22)) 

S (p[0,N-iyPN) P " = S(p+ Ar _ 1] ;p w , 7 r). □ 

Proof of Theorem 11. We also apply Proposition 17. This time we need to modify 
the state spaces: 

S 2j - 1 =Y{ b p+,p+ j _ 1] ), j = l,...,N; 
S 2k =Y( b p+ 7 p+ k _ 1] ), k = l,...,N-l. 



and 
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Also set 

W 1 = ^pt^tj-iyPi^ j = l,...,N-l; 
=P t CPo,P[ 1 , j --i];Pj), 3 = 1,...,AT- 1; 

m n — S(p^ )N _ 1 ^p N ), 

^ = Po , pjj-i] ; *)> i = i. ■ • ■ , n - 1, 

Again, the commutation relations (20) follow from Proposition 18, and the matrix 
of transition probabilities ft") from (21) coincides with The claim follows from 
(23), Proposition 17, and the relation (cf. (22)) 

Hpp,n-i]>Pn) P " = s Cpo,P[i,n-i]>Pn)- d 



10. Application to the two-sided Schur processes 

Let us put the two-sided Schur process of Definition 4 into the general framework. 
We need some notation. For n > 1, an admissible matrix M, cf. Definition 2, 
and oi, . . . , a n > in the analyticity annulus of H(M; u _1 ), define 

T A „(ai, . . . , a„; M) = l-n^ s a/ m (M), A, /x G GT„. 

i=i 

For arbitrary ai, . . . , a n > also set (A £ GT„, /U e GT„_i) 

1^/1 ixdetfaf^']":^ 
T v (a 1; . . . , a n ) = — [[ I — - - A _ ^ J *A/ M (an)- 

Thus, we have matrices T(ai, . . . , a„; M) with rows and column parameterized by 
GT„, and matrices T(a\, . . . ,a n ) with rows parameterized by GT„ and columns 
parameterized by GT„_i. 

Proposition 19. In the above assumptions, the matrices T(ai, . . . , a n ; M) and 

T(ai, . . . ,a n ) are stochastic, and the following commutation relation holds: 

T(a u ...,a n ;M) T(a u ...,a n ) = T(a u ...,a n ) T(a u . . . , o„_i; M). 

For admissible matrices Mx,M 2 and a\,...,a n > in the analyticity annuli of 
H(M i ;u^ 1 ), i = 1,2, we also have the commutation relation 

T(oi, ...,«„; Mi) T(ai, . . . ,a„; M 2 ) = T(ai, . . . ,a„; M 2 ) T(ai, . . . ,a„; Mi). 



Proof. Follows from Propositions 2.8-2.10 and Lemma 2.13(h) of [BF]. □ 
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Consider now the two-sided Schur process T(ai, . . . , ajy; M; of Definition 4. 
Set n = c(l) + • • • + c(AT) + AT, and (c(0) := 0) 

5,-=GT fc , c(fc - f) + k < j < c(k) + k, k = l,...,N. 

Define the stochastic links by 

A c(fc:l)S-i= T (oi,- ■•,«*), k = 2,...,N; 

A^:lj+J^|_ 1 =T(ai,...,a fc ;M( fc - , )) ) fc = l,...,iV, J = 1, . . . ,c(fc). 

Also define a probability distribution m* on <S„ = GTjv via 



m„(Aj- — — — Hv . A £ <url AT) 

det[>i J *j (<*<)] j>: ,- = i 

where we used the notation (15). 
These definitions imply that 

T( fll ,..,a N ;M;*)(A«,...,Am) 

= m* A «_ x ( A ( A Wc(iV-i))) _ A 2 ( A (i,i) ) A (i, )) _ 

Note that this proves formula (16) for the partition function since 

N - n — rf ^ 1 



dct 



fe=i ~ i=i 



a fe a 



Proof of Theorem 13. Once again we apply Proposition 17. We set 5, = Sj for 
j = 1, . . . , n; Aj_ 1 = Aj_j for j = 2,...,n; and m„ = m*, 

Pj- =T(ai,...,a fc ;Q*), c(fc - 1) + fc < j < c(fc) + fc, k = l,...,N. 

Note that H{Q t ]u) = H(Q:u^ 1 ) and s^/^Q*) — s^/xiQ) for signatures A and /i 
of the same length. 

The claim now follows from Proposition 17 as the needed commutativity relations 
are given in Proposition 19, and by the Cauchy-Binet identity 



det [ap^ 3 (oi)] l J=l Il"=i ff (<3; Oj) 

det \a H ~ 3 ] n 

x E ■ ; :; v , i del 'I* ,\ „ f ; : v f : r 1 ^wi^fw. 



22 



ALEXEI BORODIN 



References 

[AESW] M. Aissen, A. Edrei, I. J. Schocnberg, A. Whitney, On the Generating Functions of 
Totally Positive Sequences, Proc. Natl. Acad. Sci. USA 37 (1951), no. 5, 303-307. 

[BFP] O. Bodini, E. Fusy, and C. Pivoteau, Random Sampling of Plane Partitions, to appear 
in Combinatorics, Probability and Computing, arXiv : 0712 . 0111. 

[BF] A. Borodin and P. Ferrari, Anisotropic growth of random surfaces in 2+1 dimensions, 
arXiv: 0804. 3035. 

[BG] A. Borodin and V. Gorin, Shuffling algorithm for boxed plane partitions, Adv. Math. 
220, no. 6, 1739-1770, arXiv:0804.3071. 

[BGR] A. Borodin, V. Gorin, and E. M. Rains, q- Distributions on plane partitions, arXiv: 
0905.0679. 

[BK] A. Borodin and J. Kuan, Random surface growth with a wall and Plancherel measures 
for O(oo), to appear in Comm. Pure Appl. Math., arXiv: 0904 . 2607. 

[BO] A. Borodin and A. Okounkov, A Fredholm determinant formula for Toeplitz determi- 
nants, Integral Equations Operator Theory 37 (2000), no. 4, 386-396, math/9907165. 

[BMRT] C. Boutillier, S. Mkrtchyan, N. Reshetikhin, and P. Tingley, Random skew plane parti- 
tions with a piecewise periodic back wall, arXiv: 0912. 3968. 

[DF] P. Diaconis and J. A. Fill, Strong Stationary Times Via a New Form of Duality, Ann. 
Prob. 18 (1990), 1483-1522. 

[E] A. Edrci, On the generating functions of totally positive sequences II, J. Analyse Math. 

2 (1952), 104-109. 

[Jl] K. Johansson, Discrete polynuclear growth and determinantal processes, Comm. Math. 
Phys. 242 (2003), 277-329; arXiv: math. PR/0206208. 

[J2] K. Johansson, The Arctic circle boundary and the Airy process, Ann. Prob. 30 (2005), 
no. 1, 1-30, arXiv: math. PR/0306216. 

[K] S. V. Kcrov, Asymptotic Representation Theory of the Symmetric Group and its Appli- 

cations in Analysis, Translations of mathematical monographs, vol. 219, Amcr. Math. 
Soc, Providence, RI, 2003. 

[M] I. G. Macdonald, Symmetric functions and Hall polynomials, 2nd edition, Oxford Uni- 

versity Press, 1995. 

[Ok] A. Okounkov, Infinite wedge and measures on partitions, Selecta Math. 7 (2001), no. 1, 
57-81; math/9907127. 

[OO] A. Okounkov and G. Olshanski, Asymptotics of Jack polynomials as the number of 
variables goes to infinity, Intern. Math. Res. Notices (1998), no. 13, 641-682; arXiv: q- 
alg/9709011. 

[OP] A. Okounkov and R. Pandharipande, Gromov-Witten theory, Hurwitz theory, and com- 
pleted cycles, Annals of Math. (2) 163 (2006), 517-560; arXiv: math/0204305. 

[OR1] A. Okounkov and N. Reshetikhin, Correlation functions of Schur process with applica- 
tions to local geometry of a random 3-dimensional Young diagram, J. Amer. Math. Soc. 
16 (2003), 581-603; arXiv: math. CO/0107056. 

[OR2] A. Okounkov and N. Reshetikhin, Random skew plane partitions and the Pearcey process, 
Comm. Math. Phys. 269 (2007), no. 3, 571-609, arXiv:math. CO/0503508. 

[ORV] A. Okounkov, N. Rcshtikhin, and C. Vafa, Quantum Calabi-Yau and Classical Crys- 
tals, The Unity of Mathematics, Birkhauscr Boston, 2006, pp. 597-618; arXiv :hep- 
th/0309208. 

[Ol] G. Olshanski, The problem of harmonic analysis on the infinite dimensional unitary 

group, J. Funct. Anal. 205 (2003), no. 2, 464-524; arXiv: math/0109193. 
[PS] M. Prahofcr and H. Spohn, Scale invariance of the PNG droplet and the Airy process, 

J. Stat. Phys. 108 (2002), 1071-1106; arXiv: math. PR/0105240. 
[T] E. Thoma, Die unzerlegbaren, positive- definiten Klassenfunktionen der abzdhlbar un- 

endlichen, symmetrischen Gruppe, Math. Zeitschr. 85 (1964), 40-61. 
[VK] A. Vcrshik and S. Kcrov, Characters and factor representations of the infinite unitary 

group, Soviet Math. Doklady 26 (1982), 570-574. 
[Vo] D. Voiculcscu, Representations factorielles de type Hi de U(oo), J. Math. Pures et Appl. 

55 (1976), 1-20. 



